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Abstract. We obtain the exceptional part of the eigenspectrum of the generalised 
Rabi model, also known as the driven Rabi model, in terms of the roots of a set 
of algebraic equations. This approach provides a product form for the wavefunction 
components and allows an explicit connection with recent results obtained for the 
wavefunction in terms of truncated confluent Heun functions. Other approaches 
are also compared. For particular parameter values the exceptional part of the 
eigenspectrum consists of doubly degenerate crossing points. We give a proof for 
the number of roots of the constraint polynomials and discuss the number of crossing 
points. 


1. Introduction 

Despite it’s simplicity, the generalised Rabi model has been solved only recently 
[H [2113111115]. The Rabi model [6] describes the simplest matter-held interaction, namely 
between a two-level atom and a single-mode bosonic held. It is thus a fundamental 
textbook model in quantum optics [7|. The generalised Rabi model has hamiltonian 

H = u a)a + g ax{a) + a) + ^az + eax i (1) 

where ax and are Pauli matrices for a two-level system with level splitting A. The 
single-mode bosonic held is described by the creation and destruction operators and 
a with [a, a^] = 1 and frequency u. The interaction between the two systems is via 
the coupling g. The Rabi model has Z 2 symmetry (parity) which is broken by the 
addition of the term e ax in the generalised version of the model. This additional term 
allows tunnelling between the two atomic states. The generalised Rabi model ([I|) is also 
referred to as the driven Rabi model [8| and is relevant to the description of various 
hybrid mechanical systems mu- Although having an analytic solution, both the Rabi 
and generalised Rabi models do not appear to be integrable in general in the Yang- 
Baxter sense [10]. However, the existence of monodromy matrices in terms of Painleve 
V has now been reported [H]. 
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The generalised Rabi model ([I]) has been solved in two ways: (i) by mapping 
the problem to the Bargmann space of analytic functions DEI, and (ii) by using 
the Bogoliubov operator method [2]. Using the former approach explicit expressions 
have been obtained mu for the wavefunction in terms of confluent Heun functions 
[I^l^ Of particular relevance here is the fact that the energy spectrum of the 
generalised Rabi model, although possessing no parity symmetry, still includes both 
regular and exceptional parts. The eigenspectrum can be determined from the analytical 
solution. The exceptional parts, known as Juddian isolated exact solutions [I6], can be 
systematically found from the conditions under which the confluent Heun functions are 
terminated as hnite polynomials |1]. Our interest here is with this exceptional part of 
the eigenspectrum, which we obtain using a different approach. 

Indeed the exceptional part of the Rabi model eigenspectrum has been obtained 
using a number of different (though ultimately related) approaches]^ Each approach 
results in the simple eigenvalue expression of a shifted oscillator, however with the system 
parameters satisfying a constraint which becomes increasingly complicated for higher 
energy levels. Of most relevance here is an approach which derives a set of Bethe-like 
algebraic equations whose solutions dehne the constraint among the system parameters 

[23] . We apply this approach in Section 2 to obtain the exceptional part of the 
eigenspectrum of the generalised Rabi model ([T]), allowing an explicit connection with 
the results obtained for the wavefunction in terms of truncated confluent Heun functions 
|3]. The approach used here provides a simple product form for the wavefunction 
components in terms of the algebraic roots. We conclude by discussing the relation 
between the various approaches in Section 3. 


2. Results 


In the Bargmann realisation [25] 


X d 

a' ^ z, a —)■ — 
dz 

(2) 

the hamiltonian ([T]) reads 


H = u z-^ + ga^(z + ^ ] + Aa^ + ea^ . 
dz \ dz J 

(3) 


Following, e.g.. Ref. |1], in terms of the two-component wavefunction 

I The connection with confluent Heun functions was made earlier for the Rabi model [SKiiiin]. 

§ See, e.g., Refs [DdTKlHKlllinilllllllllllll]. A different generalised Rabi model is considered 
in Ref. |23] with both rotating and counter-rotating terms, i.e., interpolating between the Jaynes- 
Cummings and Rabi models. 
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the Schrodinger equation Hfj = Efj gives rise to a pair of coupled equations for 
and namely 


{uz + 5')^^ + { 9 z + €- E)^jJ+ + A^jJ_ = 0 , 


{uz - g) 


dfj- 

dz 


{gz + e + E)f)_ + = 0 . 


(5) 

( 6 ) 


Two sets of solutions for 'ifj^^z) and 'if-{z) can be obtained. For the hrst set, the 
substitution fj±{z) = leads to the coupled equations 


uz + g)-^ - [ — + E - e 


dz 


u 


(j)l{z) = -A(j)f{z ), 

= -A(j)\{z) 


{uz - g)-^ - (2gz - — + E + e 
dz \ u 

Eliminating 4>f{z) gives the second order differential equation 

d^fUz) 


(7) 

( 8 ) 


{uz - g){uz + g)- 


d?z 


+ 


9 


—2guz + {u — 2g — 2Eu)z + — {2g — u — 2eu) 

(jJ 

2g I ^ < TT , 1 I IT'S A 2 ,2 


+ E-e]z + E^-A^-e^ + 
u / u 


2eg^ / 


u^ 


d(j)i{z) 

dz 

i#.+(z) =0. 


(9) 


This is a case of the general second order differential equation considered by Zhang 
Applying the result of Theorem 1.1 therein [26] gives the wavefunction component 
in the factorised form 

N 

'fl{z)=e->^^/^l[{z-z,) ( 10 ) 

^=1 

where the roots Zi satisfy the set of algebraic equations (details are given in Appendix 

A) 


E 




2u 



2u‘^gzf + {2Nu — u + 2e) u'^Zi + u'^g + 2eug — 2g^ 

/ ,2^2 I 2 

^ — 9 

Nu'^ + 2eu Nu"^ — u"^ 

-+-+ 2^ 

uZi - g uZi + g 


( 11 ) 


for z = 1,..., iV. The system parameters obey the constraint 

N 

A^+ 2Ng‘^+ 2ug^Zi = t). ( 12 ) 

i=l 

The energy of these states is given by 

E = Nu- — + e. (13) 

OJ 

The corresponding wavefunction component ipfiz) can be determined using the result 
(fTOj) and equation ([7|). For e = 0 the algebraic equations (fTTD reduce to those obtained 
by the same approach [2T] . The energy expression 0131) has been given in Ref. [I] , where 
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it follows as a condition for the general solution given in terms of the confluent Heun 
functions to truncate to a polynomial with N terms. 

Another set of solutions follow from the substitution leading 

to the coupled equations 


{uz - g) 


fL 

dz 

d 


{ujz + 5') + ( ‘^gz H- E -\- e 


fl{z) = -Aff{z), 


dz 


9 


ui 


+ E + e 


if{z) = -Afl{z). 


(14) 

(15) 


Proceeding as above, these equations can be solved for the wavefunction components in 
the form 

N 

^f^z)=e^^/-l[-{z-Zi) (16) 


2=1 


where the roots Zi satisfy the algebraic equations 


2uj 


N 

V 


—2u^gzf + {2Nu — u — 2e) urzi — urg + 2eijjg + 2g^ 
lzAzI — g^ 

— 2euj 

H-:- ^9 


uzi - g ojZi + g 
for i = 1,..., iV. The system parameters now obey the constraint 

N 

+ 2Ng‘^ — 2ug ^ Zi = 0 , 


(17) 


(18) 


2=1 


with energy 

E = Nu - — - e. (19) 

OJ 

The corresponding wavefunction component 'f\{z) = e^^^^4)\{z) follows from the result 
flTB]) and equation (IT^ . The energy expression flT^ has also been given in Ref. [Ij, again 
following from the condition for truncation of the general solution given in terms of the 
confluent Heun functions. This other set of solutions was not considered for e = 0 |21j . 
The resemblance of algebraic equations of this type with Richardson BCS equations of 
Gaudin type has been noted [23] . 

There is clearly a symmetry between the two sets of solutions. Namely the algebraic 
equations flTTl) and flTTl) are equivalent under the transformation Zi -H- —Zi, e -f-)- —e. This 
corresponds to the related symmetry f)\{z,e) = 'iff^—z, —e), 'ff^z^e) = —e) in 

the wavefunction components. This symmetry is further discussed in Ref. [1] and is 
well known in the e = 0 case (see, e.g.. Ref. [I9]). The — sign has been inserted into 
equation ffT6l) to ensure this symmetry. 


2.1. Examples 

We now turn to some specihc examples. First consider N = 1. The energy is 
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and the algebraic equations flTT]) reduce to 

2u‘^gzl + (a; + 2e) u^Zi + u^g + 2eug — 2g^ = 0 . (21) 

The two solutions are 

2g^ 


Zl = 


I 

U} ' 


Ul 


2eoj 


2ug 


( 22 ) 


Substitution into the constraint flT^ gives = 0 and 

+ 2ea;, (23) 

respectively. The value A^ = 0 obtained from the hrst solution corresponds to the 
degenerate atomic limit, which we discuss further below. The second solution in 
gives the wavefunction components 

'2ujgz + + 2euj — 2g'^' 






2bjg 


'f\{z) = e 




A 


(24) 

(25) 


where in the last equation, we made use of the simplifying constraint 
On the other hand, equation flTT)) becomes 

— 200 "^gzl + (a; — 2e) oo'^Zi — u'^g + 2eujg + 2g^ = 0 . 

The two solutions are 

—2g'^ + — 2eu 


9 

Zl = 

u 


2ug 


(26) 


(27) 


Substitution into the constraint fflSj) gives A^ = 0 and A^ + d^f^ = u'^ — 2eu, respectively. 
The relevant energy and wavefunction components are 

E = u- — -e, (28) 

CJ 




^f(z) 


A 


—2ugz + — 2euj — 2g‘^ 

2ug 


(29) 

(30) 


The results for N = 1 agree with those obtained from the truncation of the confluent 
Heun functions |1], within a harmless renormalisation of the wavefunction components. 

As a further check, consider N = 2 for which equations (1TT|1 are seen to give six 
sets of solutions. For the hrst set, Zi = Z 2 = —g/oj, the constraint relation (IT^ gives 
A^ = 0. The solution zi = Z 2 = g/oj gives the unphysical constraint A^ + 8g‘^ = 0. As 
in dealing with Bethe Ansatz equations, the equations need to be solved numerically for 
hnite sizes. For the simplest case e = 0, the solutions with distinct roots have 

-5 + 4^2 ± ^9 + 8^2 + 15^4 


Zl -\- Z2 — 


^~9 


(31) 
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with g = g/uj. Substitution into the constraint relation fll2p and squaring gives the 
known result, namely + 12A^5f^ — SA^ca^ + ?>2g'^ — 32u‘^g‘^ + = 0. Equations flTTI) 

and (ITSjl give the same constraint. The explicit wavefunction components for N = 2 
and e 7^ 0 are 


{z^^aiz + al) , 
i>^_{z) = [z^ — aiz + af'") , 

where 

A2 + %2 ^ ^4 ^ g^2^2 ^ g^4 _ ^2^2 _ 2^2^^ 

“ 2gu ’ 

^ A A"* + 6A^(/^ — A^w^ — 2A^ew 

4AgV ' 


(32) 

(33) 

(34) 

(35) 

(36) 

(37) 


2.2. Degenerate atomic limit 

Some comments can be made about the degenerate atomic limit A = 0 for general N. 
The degenerate solutions Zi = —gjoj^ for z = satisfy the algebraic equations 

OB, with A^ = 0 following from the constraint relation fll2l) . The energy is given by 
flT^ with fllOl) giving the wavefunction component 

V.yj)=e-»"'“h + ^)". (38) 

\ uj / 

This is precisely the solution obtained for the equivalent displaced harmonic oscillator 
in the Bargmann space 125 ]. The related solution [18] similarly follows from equations 

dH-dH. 


3. Discussion 

It is interesting to compare the various approaches for deriving the exceptional part of 
the eigenspectrum. We have derived a set of algebraic equations (ITT]) for the exceptional 
part of the eigenspectrum of the generalised Rabi model ([I]) using a method [211126] akin 
to the functional or analytic Bethe Ansatz. Although the energies have a simple form 
fir3|l . the constraint relations flT^ and wavefunction components flTOj) are given in terms 
of the Bethe-like roots Zj. The constraint relations can be generated by a number of 
methods. It is known, for example, that the coefficients of the wavefunction components 
satisfy a system of 2N + 1 linear equations, with the constraint emerging as a condition 
for the determinant to vanish. One can also determine a recurrence relation leading to 
the constraint relations (see, e.g., [20]). 

For the generalised Rabi model considered here, in terms of the series expansion 
coefficients hn for the conffuent Heun function ^ h^x"', where x = the recurrence 
relation is 

^nhji Bjihn—l T Cnhn—2 > 


(39) 
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with initial conditions h_i = 0, ho = 1. The coefficients are given by 

An = n(— 1 + n — N — 2e/u) , (40) 

B^ = {l-n + N)^- 4(n - + 2{l-n + N)e/u: , (41) 

Cn=^{-2 + n-N)gyu\ (42) 

This result follows from Ref. [1] specified to the exceptional points[Jj] Indeed, the general 
three-term recurrence relation is central to the analytic solution of the generalised Rabi 
model [H El El SI E]. The approach taken here effectively gives the factorisation of the 
truncated confluent Heun functions at the exceptional points. 

The recurrence relation (l39|) . and in particular the condition Cn +2 = 0, ensures that 
the inhnite series expansion for the confluent Heun function terminates with hn = 0 for 
n > N. The value h^^i = 0 determines the constraint relation for given N. The hrst 


few polynomials obtained in this way are 

= 0 , (43) 

(A^ + 4^2 _ a;2 _ = 0 , (44) 

A2 + A2(12c/2 - 5uj^ - 6eu) + 32^ - 32e/a; + 8eV 

- 32u^g^ + 12ea;^ + 4a;^) = 0 , (45) 

A2(A® + 2A^(12c/ 2 - 7a;2 - Qeuj) + + 44e^u^ - 232u;^g^ 

+ 176/ + 16€u;(6u;^ - llg'^)) - 12(-32/ + 24a;/(2e + 3a;) 

- I2uj‘^g‘^{2e^ + 5ea; + 30;^) + 4e^a;^ 

+ 126 ^ 0 ;^ + lleo;^ + 3a;®) = 0 , (46) 


for A = 0,1, 2, 3, respectively. The N = 1 result fSTj) is as given in fl23|) . with fSSjl the 
example given in Ref. [1] . The constraint polynomials for given N are generated readily 
enough via the recurrence relation. A similar recurrence relation can be written down 
corresponding to the solutions = This results in the same constraint 

polynomials as given in the above examples, however with e y-)- —e. In contrast the 
approach used here gives the closed form expressions flT^ and flTSl) . albeit in terms 
of the roots of the algebraic equations (1TT|1 and m- These equations remain to be 
explored. 

3.1. Degenerate crossing points 

It has been noted that when e is an integer multiple oi uj/2 the exceptional eigenvalues 
considered here for the generalised Rabi model are crossing points in the eigenspectrum 
as a function of the coupling g na. For e = 0 these are the well known Judd points, 
for which Kus [IH] provided a proof that for each N there are N crossings in the range 

II Note that, taking oj = 1 and e = 0 for simplicity, this recurrence relation can also be written in the 
form (to -I- 1)(to — N)hm+i + [A^ -I- — [N — to)^] hm + 45^(A -I- 1 — m)hm-i = 0, which differs 

from the recurrence relations given elsewhere, e.g., in Refs [HI [20]. Presumabley this is because the 
coefficients in the recurrence relation change with the expansion variable, in this case either x defined 
above or z. 
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g 


Figure 1. The first few energy levels E in the eigenspectrum of the generalised Rabi 
model as a function of the coupling g denoted by grey (thin) lines. The parameter 
values are e = with A = 1.2 and oj = 1. For this particular value of e the exceptional 
part of the eigenspectrum consists of the doubly degenerate crossing points. The blue 
(thick) lines are the energy curves E = Nuj — /oj + e for N = 1,..., 5 



g 


Figure 2. The first few energy levels E in the eigenspectrum of the generalised Rabi 
model as a function of the coupling g denoted by grey (thin) lines. The parameter 
values are e = ^oj with A = 1.5 and oj = 1. The blue (thick) lines are the energy 
curves E = Noo — g^ /oj + e for A = 1,..., 5. 


0 < A < 1. More generally Kus established that for k < IS. < k + 1 there are N — k oi 
them. Figures [Hand [2] show the hrst few levels in the eigenspectrum of the generalised 
Rabi model for e = at different values of A. It is clear from Figure [1] that for the hrst 
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flf 


Figure 3. The first few energy levels E in the eigenspectrum of the generalised Rabi 
model as a function of the coupling g denoted by grey (thin) lines. The parameter 
values are A = 1.2 and ui = 1 with now e = 0.3. The blue (thick) lines are the energy 
curves E = Nuj — g'^/uj + e for TV = 1,..., 5. The red (thick) lines are the energy curves 
E = Nuj — g'^/oj — e for N = 1,..., 5. For this value of e there are no crossing points. 
The first nine exceptional points in the eigenspectrum are indicated by red (TV = 1), 
blue (TV = 2) and green (TV = 3) circles. 


few values of N there are N crossings for A = 1.2. We expect that this is indeed the case 
for all N in the range 0 < A < y/2. Following Kus |T8] we are able to prove a theorem 
by induction in Appendix B giving the number of roots of the constraint polynomial 
and thus the number of exceptional points for given N and e. Specihcally, dehning the 
function Qk{x) = {l/k\)Pk{x), where Pn{x) = 0 is the constraint polynomial, we have 
established the following theorem. 


Theorem. For 0 < A/u < + 2e/oj, Qk{x) has exactly k different, positive 

, (fc) (fe) (fc) 

roots a\ ,02 ,..., ; moreover 


0 < a) < a) < 02 < 02 < • • • < < o^ 


(47) 


where denote the roots of Qk-i{x). 


More generally, as discussed in Appendix B, one can prove that there are N — k 
roots of the constraint polynomial for A in the range 

+ 2ke/oj < A/oj < \/{k + 1)^ + 2{k + l)e/a;. (48) 


Indeed, in Figure |2] we see that for the hrst few values of N there are N — 1 crossings 
for A = 1.5 when e = |a;. However, to prove the number of level crossings for 
e = requires a further step. Figure [3] shows the hrst few exceptional points in 
the eigenspectrum at the typical value e = 0.3 a; where there are no crossing points. It 
was pointed out how these exceptional points merge to form doubly degenerate crossing 
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points as e —)■ |ci; [1]. For the number of level crossings, it is necessary to prove that the 
two corresponding sets of roots of the constraint polynomials coincide for given N when 
e = see Appendix B. This does not seem so straightforward to prove however, for 
general N. Nevertheless, we are confident that there are indeed N — k doubly degenerate 
crossing points in the range fHHl) . It would be fascinating, though seemingly unlikely, if 
such theorems could be proved via the algebraic equations obtained in this paper. 
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Appendix A. Derivation of the algebraic eqnations 

The second order differential equation ([9]) is of the general form 

R(z) = 0, 


, , df , . d 


where 


X{z) = J2<^kX^ Y{z) = J2hz^ Z{z) = J2< 


= 2_^CkZ 

k=0 k=0 k=0 

Comparing with equation ([9]), the nonzero coefficients are 
ao = - 5'^ a2 = uj‘^ , 

2(72 

bo = - 2eg — tzg , bi = — 2g^ — 2Ee , 62 = —‘^<xg , 

u 

r -,2 A 9 9 g"^ 2g^ ^ 

cq = E — A — e H- - , C 2 =-h 2Eg — 2eg . 


OJ OJ^ 


00 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 


Zhang’s Theorem 1.1 [2^ states that (lA.ip has a degree n polynomial solution 

n 

*^(^)= 


2=1 


with distinct roots zi,Z 2 , ■ ■ ■, Zn- The values of the coefficients cq, ci, C 2 are given by 

C 2 = — n{n — 1)04 — nbo , (A-7) 


Cl = - [ 2 {n - 1)04 + 63] ^ Zi - n{n - 1)03 - nb2 , 

2=1 

n n 

Co = - [ 2 {n - 1)04 + 63] ^ Zi - 204 ZiZj 

2=1 i<j 

n 

- [ 2 {n - 1)03 + 62] ^ Zi - n{n — 1)02 - nbi. 


(A.8) 


(A.9) 


2=1 
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The roots zi, Z 2 , ■ ■ ■, Zn satisfy the algebraic equations 
2 , + b 2 zf + biZi + 60 


n 




Zi - Zj a^zf + aszf + a2zf + aiZi + oq 


= 0, 


(A. 10) 


j¥^i 

ioY i = 1,2 ,... ,n. 

Substituting the values flA.3P - flA.5j) into flA.7D - flA.10j) gives 0 = 0, the energy 
expression m, the constraint flT^ and the algebraic equations flTT]) . respectively. In a 
similar fashion we arrive at equations fITTD . flT8|) and ([T9]). 


Appendix B. Proof for the number of exceptional points 


To prove the number of roots of the constraint polynomial it is convenient to generalise 
the recurrence relation obtained by Kus [TS] rather than the recurrence relation fl39p . 
For nonzero e we use the recurrence relation 

Pq = 1 , Pi = 4:g‘^ + — 2eu , 

Pk = [k{2gY -I- — 2keuj] Pk-i 

- k{k - l){n - k + l){2gfu‘^Pk-2 ■ (B.l) 

The equation Pk = 0 when k = N dehnes the constraint polynomial, which can be 
written here in the form 


[N{2gf + - AV - 2Neoj\ P^-i - N{N - l)(2c/) VPjv -2 = 0.(B.2) 

We now £x the value of N and set x = {2gY. Thus 
Qo{x) = l, Qi{x)=x-ax, 

Qk{x) = (x- ak)Qk-i(x) - (dkX Qk-2{x) , (B.3) 


where Qk(x) = (l/k!)Pk(x) and 

ak = (k^oo^ + 2k€0L> — A^)/k, (B.4) 

/3k = (n - k + l)uj^ . (B.5) 


Now, following Kus [18], we can prove the theorem stated in section 3.1. 

Proof. For 0 < A/cu < -^1 -|- 2e/a; we have > 0 with always /3k > 0. From the 
dehnitions we also have = ai > 0 and Q2{x) = {x — 0)2)(x — — / 32 X. Thus 

<52(0) = Q!2ai^^ > 0, Q2{a^i^) = -/ 32 a^i^ < 0 and sgn (52(00) = 1 , where 

f -1 a < 0 

0 a = 0 . 

1 a > 0 


These results and relations flTTD prove that Q 2 {x) = {x — a\ ){x — 02 ) and 0 < a) < 

The general proof now proceeds by induction. Assume that the theorem is valid 
for i < k, i.e.. 


Qk-i{x) = {x- af ^^)---(x- , 

Qk- 2 (x) = (x- af • • • (x - , (B.6) 
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and 


,(^-1) / .(fc-2) / Ak-1) (k-2) 

Cio 


0 < a\ ’ < a\ ’ < 02 
Then from flB.3D we have 

Qk{x) = {x-ak){x- ■■■{x- 4 -”/^) 

- f3kx{x - of ■ ■ ■ (x - a^k- 2 ^) ■ 


^ ^ {k-2) ^ (fe-1) 

< ■ ■ ■ < ^k -2 < ®fc-l 


Thus 


(B.7) 


(B.8) 


sgnQfc(O) = sgn(-l)^afcaf ^^4^" = {-if , 

sgnQfc(af“^^) = - sgn af~^^ {af~^^ - - afsf^] 

^)) = {-lf-\ 


x(af - of 
sgnQfc(oo) =1. 


(fc-l) {k-2) 


(B.9) 


This implies that 

Qk{x) = {x - ■■■(x - , (B.IO) 

and ..., fulhll the inequalities flT7|l . The theorem is thus proved. 

Again following Kus [H] similar theorems can be proved for other values of A. The 
ranges of A follow from the values of k for which ak dehned in fIB.dp is no longer positive. 
In this way one can prove that there are N — k roots of the constraint relation in the 
range fH8l) . When e = |a; there are thus expected to be A —A; roots (exceptional points) 
in the range \/k? + k < A/u < {k + If + k + 1 . There are still expected to be A — fc 
points for given A at A = fkP' + k. However, for this special value we observe that at 
the left most edge of the energy plots two energy levels merge, rather than crossing, so 
technically they are not crossing points. 

To prove the double degeneracy at the crossing points we need to consider the other 
set of recurrence relations for nonzero e, which can be written as 


Fq = 1, F( = Ag'^ + A^ — — 2e a;, 

F^ = [k{2gf + A^ - A;V + 2keu\ F^^ 

-k{k-l){n-k + 1)(2(7) VF ^2 • (B.ll) 


The equation F^ = 0 when k = N dehnes the constraint polynomial. We can repeat 
the above working to arrive at similar results for the number of roots and thus number 
of exceptional points. Now one can prove that the function Qf{x) = (l/fc!)F^(a;) has 
N — k different positive roots for A in the range 

a/A;^ — 2 kejuj < A/oj < \/{k + 1)^ — 2{k + l)e/a;. (B.12) 

When A: = 0, Qf{x) has A roots for 0 < A/o; < \/l — 2e/a;. We thus have that for 
< A/uj < \/{k + 1)^ — 2 {k + l)e/a; both Qn{x) and Qf{x) have N — k 
roots. In particular, for \/{k + 1)^ — 2{k + l)e/a; < A/uj < \/{k + 1)^ + 2{k + l)e/a;, 
Qn{x) has N — k roots and Qf{x) has A —A; —1 roots. Precisely at e = the intervals 
+ ‘^ke/u, \/{k + 1)2 - 2 {k + l)e/ vanish. This implies that for arbitrary A, 
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Qn{x) has N — k roots and has N — k — 1 roots. This is the situation seen, 

for example, in Figure 3. We have been able to show numerically that the roots of 
the polynomials Qn{x) and coincide for all values of N when e = The 

exceptional points are thus doubly degenerate at the crossing points. However, we have 
not so far been able to prove this analytically. 
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